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Theoretical  Studies  of  Kinetic  Mechanisms 
of  Negative  Ion  Formation  in  Plasmas 


ABSTRACT 

This  technical  program  constitutes  a  theoretical  research  investigation  of 
the  kinetic  mechanisms  of  negative  ion  formation  in  plasmas.  This  study  was 
directed  toward  elucidating  the  mechanisms  of  the  most  important  volume  dependent 
reactions  that  occur  in  hydrogen-ion,  H~(D"),  source  devices,  primarily  of  the 
Belchenko-Dimov-Dudnikov  (BDD)  type  and  toward  evaluating  other  light  negative 
anions,  6uch  as  Li",  as  possible  sources.  The  primary  goal  of  this  research 
program  was  to  identify  the  most  important  reactions  leading  to  negative  ion 
production  or  destruction  and  to  estimate  these  reaction  rates  as  a  function  of 
system  parameters  such  as  density,  composition  and  temperature.  A  further  goal 
was  to  explore  new  chemical  sources  for  the  production  of  light  mass  negative 
atomic  ions.  The  results  of  thi6  program  furnish  data  and  provide  direction  for 
more  detailed  investigations  into  the  kinetics  of  both  gas  phase  and  gas-surface 
reaction  rates  of  importance  in  ion  source  devices  and  provide  input  for  reliable 
modeling  of  such  systems.  This  investigation  was  carried  out  using  quantum 
mechanical  methods.  Both  ab  initio  and  density  functional  approaches  were 
employed  in  these  studies. 
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STATEMENT  OF  WORK 


The  contractor  shall  conduct  a  theoretical  research  investigation  of  the 
electronic  structure  of  negative  molecular  ions  and  the  kinetic  mechanisms  of 
negative  ion  formation  in  plasmas.  This  investigation  shall  include  the 
calculation  of  potential  energy  hypersurfaces  for  ion-atom  and  ion-molecule  gas 
phase  interactions  that  occur  in  magnetron-type  negative  ion  source  devices.  The 
contractor  will  perform  theoretical  studies  of  H  surface  adsorption  to  reduce  the 
surface  conversion  mechanism  for  H~  production  into  a  series  of  well 
characterized  elementary  processes.  In  addition,  calculations  shall  be  performed 
to  estimate  the  rates  of  the  most  important  volume  dependent  reactions  leading  to 
H  production  or  destruction  in  such  devices.  Experimental  data  collected  at 
LBL,  LASL  and  AVCO-Everett  from  their  ongoing  studies  with  negative  ion  source 
devices  shall  be  iteratively  examined  and  integrated  with  the  results  generated 
from  this  theoretical  investigation.  The  computations  shall  be  carried  out  with 
ab  initio  and  density  functional  methods  using  electronic  structure  and  kinetic 
cross-section  codes  that  have  previously  been  developed  or  modified  at  this 
Center.  Specifically,  this  investigation  shall  comprise  the  following: 

a.  Quantum  mechanical  calculations  of  the  kinetic  cross-sections  for  volume 
dependent  reactions  leading  to  the  collisional  formation  and  loss  (destruction) 
of  H". 


b.  Calculation  of  dissociative  attachment  cross  sections  for  e  ♦  Li2 

■*  Li  *  Li  as  a  function  of  collisional  energy  and  rovibrational  energy  of  the 
Li2  molecule. 

c.  Development  of  a  computer  code  and  calculation  of  the  electronic 
excitation  cross  sections  for  Li2  via  electron  collisions. 

d.  Continued  studies  of  negative  ion  formation  on  surfaces  using  large 
cluster  representations  with  fractional  hydrogen  coverage. 


I  »■  Hj  1 1  m  1 1  9 
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BACKGROUND 


Recent  successes  in  the  Tokamak  program  and  in  other  controlled  thermo¬ 
nuclear  research  programs  (Ref.  1)  have  focused  attention  on  the  problem  of 
developing  an  efficient  high-energy  particle  beam  source.  For  several  applica¬ 
tions,  a  neutral  beam  with  energies  above  200  keV  is  desired.  The  acceleration 
of  negative  ions  (primarily  H-  or  D-)  to  such  energies,  followed  by  neutraliza¬ 
tion  through  a  stripping  reaction,  appears  at  the  present  time  to  be  the  most 
efficient  approach  for  producing  a  high  energy  neutral  particle  beam. 


In  another  area,  the  feasibility  of  particle  beams  (both  charged  and 
neutral)  as  military  weapons  has  been  under  study.  The  proposed  endoatmospheric 
applications  require  high  intensity  sources  and  mainly  focus  on  the  problems  of 
beam  stability  and  propagation  characteristics.  Exoatmospheric  applications 
require  lower  intensity  sources  but  of  very  high  quality.  Design  goals  are 
highly  collimated  beams  with  a  narrow  energy  spread. 

A  magnetron-type  negative  H“  source  has  been  reported  by  Belchenko,  Dimov 
and  Dudnikov  (BDD)  (Ref.  2)  that  has  produced  H“  current  densities  of  several  A 
cm-2.  This  device  operates  as  a  plasma  discharge  in  an  atmosphere  of  hydrogen 
gas  with  cesium  or  other  alkalis  present  at  ~  0.01  percent.  The  mechanism  for 
the  production  of  H“  is  believed  to  involve  a  catalyzed  surface  reaction  whereby 
H-  ions  are  produced  by  transfer  of  electrons  from  Cs  atoms  that  are  adsorbed  on 
the  cathode  surface  (Refs.  3,  4).  The  detailed  kinetic  mechanisms  of  such 
surface  reactions  are  still  uncertain  and  parametric  experimental  studies  are 
currently  underway  at  LASL,  LBL,  IRT  and  at  other  laboratories  to  elucidate  the 
mechanisms  and  operating  characteristics  of  BDD  and  similar  devices  and  to 
develop  information  for  their  scale-up  to  higher  current  densities.  Diagnostics 
of  H-  source  devices  are  also  underway  at  Brookhaven  (Ref.  5)  using  beam  probe 
I  and  spectroscopic  techniques. 

Concurrent  with  these  surface-plasma  reactions  are  several  electron-H 
>  volume-dependent  processes  that  may  lead  either  to  the  production  of  H”  or,  in 

reverse,  may  act  as  important  destructive  processes  of  the  negative  ions  after 
they  are  formed.  Photodetachment  experiments  by  Bacal  and  Hamilton  (Ref.  6)  in 
hydrogen  plasmas  indicate  H-  densities  100  times  larger  than  that  predicted  from 
;  simple  electron  attachment  mechanisms.  Further,  these  experiments  indicate  a 

nonlinear  dependence  of  the  production  of  H”  on  electron  density,  at  least  for 
i  densities  less  than  1010  cm-3.  Several  mechanisms  (Ref.  7)  have  been  proposed  to 

'  explain  these  volume-dependent  H“  production  processes  but  the  cross-sections  for 

many  of  the  reactions  leading  to  H~  are  uncertain. 

By  a  large  margin,  the  major  thrust  of  the  development  of  particle  beam 
technology  has  been  experimental  in  nature.  Theoretical  studies  of  surface 
!  catalyzed  reactions  leading  to  H“  production  have  been  undertaken  by  LLL 

(Ref.  4),  but  only  a  relatively  small  theoretical  effort  (Ref.  8)  has  been 
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directed  toward  detailed  calculations  of  cross-sections  for  volume-dependent 


reactions.  Because  of  the  difficulty  of  conducting  experiments  to  measure  these 
unknown  cross-sections,  the  development  of  a  sound  theoretical  method  for 
calculating  these  state-to-state  processes  appeared  desirable.  Although 
relatively  little  work  of  this  nature  has  been  done  in  the  past,  enough 
theoretical  work  is  available  to  indicate  that  the  development  of  such  procedures 
can  be  made  practical,  particularly  if  good  wave  functions  and  potential  energy 
surfaces  are  available  for  the  interacting  species. 

The  original  goal  of  this  technical  program  was  to  delineate  the  most 
important  volume  dependent  reactions  that  lead  to  H"(D”)  production  (or 
destruction)  and  to  calculate  the  energetics  and  cross-sections  of  such  reactions 
as  a  function  of  system  parameters.  Attention  was  also  directed  toward  those 
kinetic  processes  that  are  of  importance  in  high  pressure  ion  source  devices, 
particularly  of  the  BDD  type.  It  was  intended  that  this  theoretical  program 
should  compliment  experimental  studies  with  negative  ion  source  devices.  A 
cooperative  program  of  ongoing  experimental  research  and  our  theoretical  studies 
are  of  mutual  benefit  and  result  in  a  better  understanding  of  the  fundamental 
physical  processes  occurring  in  ion  source  devices. 

A  further  goal  of  this  research  program  is  to  explore  new  chemical  sources 
for  the  production  of  light  negative  atomic  ions,  in  particular  H”  and  Li". 
Current  studies  include  the  formation  of  H  or  D-  via  photodissociation  of 
unsaturated  hydrocarbon  negative  ions  and  the  formation  of  Li-  via  dissociative- 
attachment  of  electrons  to  molecular  Li2.  This  latter  process  appears  feasible 
based  on  our  current  knowledge  of  the  energetics  of  this  system.  Our  recent 
research  efforts  have  been  directed  toward  three  general  areas:  mutual 
neutralization  and  collisional  dissociation  reactions  in  negative-ion  plasmas, 
analysis  of  dissociative-attachment  in  the  e  +  Li2  system,  including  the 
calculation  of  electron  excitation  cross  sections  for  e  +  Li2,  and  further 
studies  of  surface  catalyzed  reactions  for  the  production  of  H-.  The  negative- 
ion  plasma  studies  entail  analysis  of  the  kinetics  of  volume  reactions  leading  to 
the  destruction  of,  for  example,  H  and  Li  .  The  e  ♦  Li2  dissociative  attachment 
studies  parallel  the  work  of  Bardsley  and  Wadehra  (Refs.  8,9)  on  the  e  +  H2 
reaction  and  are  being  carried  out  in  collaboration  with  Professor  Wadehra. 
Although  a  phenominological  model  for  negative-ion  formation  on  metal  surfaces 
has  been  proposed  by  Hiskes,  et  al.  (Ref.  4),  a  basic  understanding  of  the 
surface  state  chemistry  has  yet  to  be  developed.  Our  work  in  this  area  includes 
density  functional  calculations  of  tungsten  and  molybdenum  clusters,  analyzed  as 
base  metals  and  with  the  addition  of  cesium  and  hydrogen  on  the  surface. 
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During  the  past  contract  year,  we  have  undertaken  studies  of  several  light 
negative  ion  systems.  As  a  continuation  of  our  work  on  Li-  formation,  we  have 
suggested  a  new  kinetic  route  for  the  formation  of  quiescent  Li-  through  Rydberg 
atom  collisions.  This  new  mechanism  is  now  under  study  by  M.  McGeoch  at  AVCO  as 
an  alternative  to  dissociative  attachment  of  e  +  Li2  as  a  source  for  Li-.  In 
addition  to  our  Li-  studies  we  have  continued  an  analysis  of  the  kinetics  of  Hej 
autodetachment.  This  work  is  being  done  in  collaboration  with  J.  Peterson  and 
Y.  Bae  at  SRI  and  has  resulted  in  several  publications  and  joint  technical 
papers.  Finally  we  have  examined  the  negative  hydrogen  ion  cluster,  H^,  to 
determine  whether  this  system  is  thermodynamically  stable  relative  to  H“  +  H2 .  I 
summary  of  our  research  is  given  below.  Four  technical  papers  and  seven 
presentations  have  resulted  from  the  past  year's  research  efforts. 


Li' 


The  negative  ions  of  light  atoms  are  currently  being  studied  for  their 
possible  application  in  gaseous  discharges,  fusion  plasmas  and  gas  lasers  (Ref. 
10).  One  source  for  the  volume  production  of  atomic  anions  is  the  process  of 
dissociative  electron  attachment  to  molecules  (Ref.  11).  This  process  has  been 
studied  experimentally  for  low  energy  electron-hydrogen  molecule  collisions  by 
Allen  and  Wong  (Ref.  12).  A  parallel  theoretical  study  has  been  reported  by 
Wadehra  and  Bardsley  (Ref.  13).  More  recently,  McGeoch  and  Schlier  (Ref.  14) 
have  examined  dissociative  attachment  (DA)  in  electron-lithium  molecule 
collisions  and  have  found  large  DA  rates  for  attachment  to  highly  vibrat ionally 
excited  Li2  molecules.  The  effect  of  vibrational  excitation  on  the  DA  rates  has 
been  studied  by  Hiskes  (Ref.  15)  for  e  +  H2  collisions  and  by  Wadehra  and 
Michels  (Ref.  16)  for  the  e  +  Li2  system. 


In  the  e  +  Li2  system,  it  has  been  shown  (Ref.  16)  that  excitation  of  Li2  to 
the  A  *  I*  state,  either  by  electron  impact  or  by  photon  pumping,  results  in  an 
enhanced  vibrational  distribution  upon  radiative  decay  of  the  A  *£*  state  to  the 


ground  X  *Eg  state  of  Li2 .  The  cross  sections  for  vibrational  excitation  of  Li2 
X  via  low  energy  electron  collisional  excitation  of  the  A  *1*  state  exhibit  a 


relatively  flat  vibrational  distribution  in  the  range  3  v"  9.  This  suggests 
that  Li2  molecules  which  are  vibrat ionally  excited  to  the  DA  threshold  for  Li" 
formation  should  exhibit  the  largest  attachment  cross  sections,  a  result  that  is 
indicated  by  the  studies  of  McGeoch  and  Schlier  (Ref.  14). 


An  alternative  mechanism  for  negative  ion  formation  was  suggested  by  Lee  and 
Mahan  (Ref.  17),  who  proposed  that  ion-pair  formation: 


X*  +  X 


X+ 


(1) 


should  occur  for  systems  where  the  energetics  for  ion-pair  formation  are 
competitive  with  associative  ionization: 
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X*  ♦  X  ♦  Xj  ♦  e.  (2) 

The  alkali  metals,  including  lithium,  are  therefore  good  candidates  for  such  a 
study  of  ion-pair  formation.  Ciocca,  et.  al.  (Ref.  18)  have  recently  reported  on 
the  formation  of  Na~  ions  in  highly  excited  (n“7-40)  Rydberg  atom  collisions. 

Our  studies  were  directed  toward  examining  the  possible  mechanisms  for  ion-pair 
formation  in  lithium  atom  collisions. 

In  order  to  determine  the  possible  paths  for  Li  excited  state  interactions, 
a  series  of  ab  initio  calculations  of  the  electronic  states  of  Li2  up  to  the 
Li[2S(3s)]  +  Li[2S(3s)]  dissociation  limit  was  undertaken.  Since  the  formation 
of  the  ion-pair  Li"*" [  2S ]  +  Li- [ 1 S ]  can  occur  only  for  *Eg  and  symmetries,  our 
studies  were  restricted  to  these  two  representations. 

Theoretical  potential  curves  for  the  *Eg  and  lj*  symmetries  of  Li2  were 
obtained  from  valence  configuration  interaction  (VCI)  calculations.  A  64 
function  Slater  orbital  basis,  containing  a  (2slp)  optimized  valence  basis 
augmented  by  a  (3s3p3d2f)  Rydberg  basis,  was  used  in  this  study.  The  basis  was 
transformed  to  symmetry  orbitals  and  Cl  calculations  were  performed  in  the 
space  of  all  symmetry  adapted  configurations  having  the  loa  orbital  doubly 
occupied.  The  resulting  Cl  expansions  contained  344  and  331  configurations  for 
the  JIg  and  ^Eu  symmetries,  respectively.  As  the  interactions  in  this  system  are 
very  long-ranged,  calculations  were  performed  over  a  range  of  internuclear 
separations  from  5  to  50  bohrs.  The  results  of  our  theoretical  calculations  are 
shown  in  Fig.  1  and  Fig.  2  for  the  and  symmetries,  respectively.  The 
ionic  curves  which  dissociate  to  Li+[*S]  ♦  Li"PS(2s2)]  are  clearly  evident  in 
both  Fig.  1  and  Fig.  2,  where  they  exhibit  a  nearly  diabatic  crossing  behavior 
with  the  normal  valence  states  of  Li2  at  large  internuclear  separations,  but  show 
strong  mixing  at  distances  <  10A.  Ion-pair  formation  may  occur  for  long-range 
interactions  of  the  type 

Li  [2S(2s)  ]  ♦  Li**  (2S(n=5)  ]  -  Li+  US]  +  Li"HS]  (3) 

but  clearly  there  are  few  effective  curve  crossing  channels  leading  to  Li+  + 

Li".  The  curve  crossing  mechanism  suggested  by  Ciocca,  et  al.  (Ref.  18),  thus 
appears  not  to  be  an  important  mechanism  for  ion-pair  production  in  the  Li2 
system. 

A  second  mechanism  suggested  by  Ciocca  et  al.  (Ref.  18)  for  alkali  negative 
ion  formation  would  involve  excited  state  negative  ion  formation  via 

Li*l2P(2p)  ]  ♦  Li**(n>9)  ♦  Li+tlS]  +  Li" [  lp( 2s2p)  ]  (4) 

or  from  electron  attachment  to  a  Rydberg  state  of  Li  followed  by  radiative 
stabilization  to  the  ground  Li~(1S(2s2)]  state.  Both  of  these  processes  can  be 
ruled  out  for  the  Lij  system  since  the  autoionization  lifetime  of  the  Li" 
[1P(2s2p)]  state  is  known  to  be  very  short  (Ref.  19)  relative  to  radiative 
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stabilization  and  higher  autoionizing  states  should  exhibit  even  shorter 
1 ifet iraes . 


The  analysis  of  our  calculated  results  given  above  leaves  us  with  radiative 
stabilization  of  excited  state  Li2  interactions  as  the  most  probable  mechanism 
for  ion-pair  formation  in  this  system.  Three  separate  types  of  interactions 
leading  to  ion-pair  production  can  be  identified: 


Li[2S(2s)]  +  Li  (n>5)  +  Lij  -*■  Li+  +  Li-  +  hv 


(18. Op  >  X  >  2. Op) 


Li*[2P(2p)]  +  Li**(n>3)  ♦  Li.**  +  Li+  ♦  Li-  ♦  hv 


(2.8p  >  X  >  0.5y) 


Li**(n>3)  ♦  Li**(n>3)  ♦  Li**  ♦  Li+  ♦  Li“  ♦  hv 


( X  <  0. 6p) 


The  mechanism  illustrated  by  Eq.  (7)  exhibits  a  quadratic  dependence  on  the 
Rydberg  atom  concentration  and,  similar  to  the  studies  reported  by  Ciocca,  et 
al .  (Ref.  18)  for  Na,  should  be  distinguishable  from  the  reactions  given  by  Eq. 
(5)  and  (6).  The  latter  reactions  should  exhibit  a  linear  dependence  on 
reaction  rate  with  Rydberg  atom  concentration.  A  preliminary  analysis  of  the 
radiative  decay  of  these  high-lying  Lij  states  shows  a  strong  preferential 
coupling  to  the  Li+[1S]  ♦  Li“[1S(2s2)]  dissociating  channel.  This  is  as 
expected  since  the  transition  moments  are  very  large  for  covalent  Rydberg 
interactions  radiating  to  a  lower  ionic  molecular  state.  In  addition,  the  range 
for  these  Rydberg  atom  interactions  is  predicted  to  be  very  large  (>30A)  where 
the  radiative  moments  are  increasing  proportional  to  the  internuclear 
separation. 


An  experiment  similar  to  that  reported  by  Ciocca  et  al  (Ref.  18)  has  been 
suggested  by  M.  McGeoch  at  AVCO  to  test  this  radiative  mechanism  for  Li+  -  Li~ 
ion-pair  formation  through  Rydberg  atom  collisions.  The  two  laser  excitation 
mechanism: 


Li[  2S(2s)  ]  +  hv  -*■  Li  l2P(2p)]  A  -  671.0  nm 


Li  ( 2P( 2p) ]  ♦  hv'  +  Li**  l 2S( 3s ) ]  A’  -  812.8  nm 


should  produce  Rydberg  Li  atoms  in  either  the  2S(3s)  or  2D(3d)  states,  provided 
near-6aturat ion  is  achieved  in  the  first  resonant  line  excitation.  The 
subsequent  reactions  of  these  laser-produced  Rydberg  state  Li  atoms  is 
illustrated  in  Table  I.  Two  photon  excitation  by  the  first  laser,  and  excitation 
of  the  2S(3s)  Rydberg  state  by  the  671  nm  laser  source,  are  both  below  the 
ionization  threshold  for  Li,  thus  eliminating  any  first  order  ionization  process. 
The  suggested  experiment  will  look  for  Li+  -  Li“  ion-pair  formation  accompanied 
by  excimer-like  radiation  at  wavelengths  <  600  nm.  The  ion  pairs  could  be 


WW 
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electrostatically  extracted  from  the  laser  reaction  region  and  mass  analyzed  to 
identify  the  Li"  signal. 


As  mentioned  above,  the  reaction  given  by  Eq .  (7)  should  exhibit  a  quadratic 
dependence  on  the  Rydberg  state  density  and  a  very  large  cross  section  for  ion- 
pair  formation.  Similar  volume  processes  to  that  represented  by  Eqs .  (5)-(7)  are 
energetically  possible  for  H*  -  H”  ion  pair  formation,  provided  that  a  useful 
mechanism  for  n  ■  2  production  in  H  can  be  found.  Very  recently  R.  McFarlane 


(Hughes)  has  developed  a  Lyman  alpha  (Lq)  laser  source  which  can  be  used  to  pump 


the  n  *  2  transition  in  H.  An  analysis  and  possible  collaboration  with 
R.  McFarlane  in  studying  this  process  is  in  progress. 


He: 


In  recent  studies  of  the  newly  discovered  (Ref.  20)  metastable  He2  anion, 
autodetachraent  to  the  He,,  [X1!*]  ♦  e  repulsive  continuum  was  observed  (Refs.  20- 


2  i  gJ  r .  . . 

22).  A  preliminary  analysis  (Ref.  21  and  22)  of  the  measured  intensity 

distribution  of  this  autodetachment  process  yields  a  vertical  separation  of 
15.7  ±  0.2  eV  between  the  metastable  He2  [4IIg]  state  and  the  repulsive  ground 
[X1!*]  state  of  He..  This  energy  separation  is  -0.5  eV  larger  than  that 


[ X 1 Zg  ]  state  of  He2-  This  energy  separation  is  -0.5  eV  largei 
calculated  using  our  best  estimate  (Ref.  23)  of  the  spectroscopic  constants  for 


Hej,  and  the  ground  state  potential  for  He2  as  represented  by  the  recent  Monte 
Carlo  calculations  by  Ceperley  and  Partridge  (Ref.  24). 


In  order  to  determine  the  source  of  this  discrepancy,  we  have  analyzed  the 
intensity  distribution  of  the  Hopfield  helium  continuian  reported  by  Huffman,  et 
al.  (Ref.  25).  This  continuum  emission  arises  from  the  transitions, 

He2  [A1!*]  ♦  [x1!*]  and  He*  [°lj:u]  *  [  X 1  Zg  ] ,  where  several  vibrational  levels  of 
the  radiating  states  contribute  owing  to  the  nature  of  the  discharge  formation  of 
the  He2  excimer.  The  spectroscopic  data  for  both  the  AlE*  and  D1!*  states  of  He2 
are  accurately  known  from  the  work  of  Tanaka  and  Yoshino  (Ref.  26)  and  Ginter 
(Ref.  27),  leaving  any  uncertainty  in  the  analysis  of  the  Hopfield  continuum 
solely  in  the  repulsive  character  of  the  ground  state  potential  of  He2 .  Since 
the  equilibrium  separation  of  the  A1!*  state  of  He2  (1.0406  A)  is  close  to  the 
predicted  value  for  the  4IIg  state  of  He2  (1.064  A),  an  analysis  of  the  He2 
continuum  intensity  should  be  a  sensitive  test  of  the  ground  state  potential  in  a 
region  where  discrepancy  between  theory  and  experiment  has  existed  (Refs.  28-30). 


The  problem  of  calculating  the  bound-free  emission  spectra  for  a  diatomic 
molecule  has  been  previously  treated  in  detail  (Refs.  31  and  32).  For  an 


analysis  of  the  spectral  intensity  distribution  arising  from  the  *  X1!*  and 


♦  X1!*  ft*0* it  ions  in  He2 ,  accurate  numerical  wavefunct ions  are  required  for 
both  the  bound  vibronic  levels  of  the  upper  radiating  states  and  the  continuum 
dissociating  state.  Such  wavefunct ions  describing  the  nuclear  motion  of  the 
system  can  be  accurately  generated,  provided  reliable  potential  energy  functions 
are  available  for  the  electronic  states  involved  in  the  transition. 
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For  the  bound  A1!*  and  D1!*  states,  accurate  RKR  potential  energy  functions 
were  generated  based  on  the  latest  analysis  of  the  spectroscopic  constants  given 
by  Huber  and  Herzberg  (Ref.  33).  For  the  ground  state  potential  of  He2»  we  have 
used  three  different  functional  forms.  The  first  was  the  exponential  fit  derived 
by  Foreman,  et  al.  (Ref.  34),  based  on  the  their  careful  analysis  of  total 
scattering  cross  sections  for  He-He  collisions.  The  second  potential  employed 
was  the  ab  initio  result  recently  reported  by  Ceperley  and  Partridge  (Ref.  24), 
based  on  quantum  Monte  Carlo  calculations  at  short  internuclear  separations 
(0.5  A  <  R  <  1.8  A).  This  theoretically  derived  ground  state  potential  is 
reported  to  be  accurate  to  ±  0.05  eV  at  0.5  A  and  ±  0.01  eV  at  1.8  A.  The  third 
potential  is  that  derived  by  Leonas  and  Sermyagin  (Ref.  35)  from  differential 
elastic  scattering  cross  section  measurements.  For  internuclear  separations 
greater  than  1.8  A,  which  affect  the  normalization  of  the  continuum  nuclear 
wavefunct ions ,  we  used  an  accurate  ab  initio  potential  that  has  been  reported  by 
Liu  and  McLean  (Ref.  36).  This  latter  potential  was  smoothly  joined  to  the 
Foreman,  Ceperley  or  Leonas  potentials  in  the  region  of  1.6  -  1.8  A. 

We  have  evaluated  the  Einstein  A-coeff icients  for  Hej  corresponding  to  the 
A1!^  ♦  X1!*  *nd  D1!*  X1!*  transitions  for  the  wavelength  region  (60  nm  <  A  < 

110  nm) ,  which  covers  the  experimental  He2  continuum  emission  reported  by 
Huffman,  et  al .  (Ref.  25).  Using  the  A-coeff icients  as  a  direct  measure  of  the 
observed  continuum  intensity,  we  show  in  Fig.  3  a  comparison  of  the  data  of 
Huffman,  et  al .  (Ref.  25)  with  our  calculated  intensity  profiles.  Our  calculated 
profiles  are  the  result  of  a  least-squares  fit  of  the  A1!*  ♦  X1!*  and  D1!*  ♦  X1  E* 
A-coeff ic ients  assuming  contributions  from  the  first  four  vibrational  levels  of 
both  transitions. 

The  calculated  vibrational  population  of  the  A1!*  and  D1!*  states  are  given 
in  Table  II  for  the  three  different  ground  state  He2  potentials  that  were 
studied.  Overall,  both  the  Foreman  (Ref.  34)  and  Ceperley  (Ref.  24)  potentials 
exhibit  a  good  fit  to  the  data  of  Huffman,  et  al.  (Ref.  25).  The  potential  given 
by  Foreman,  et  al .  (Ref.  34)  yields  a  statistically  better  fit  with  *2  ■  0.002  as 
compared  to  y2  *  0.01  for  the  Ceperley  and  Partridge  potential  (Ref.  24).  A 
significantly  poorer  fit  with  y2  ■  0.47  is  found  for  the  potential  proposed  by 
Leonas,  et  al.  (Ref.  35).  All  other  ground  state  potentials  that  we  have 
explored,  including  those  reported  by  Jordan  and  Amdur  (Ref.  37),  yield 
significantly  poorer  statistical  fits.  We  have  not  analyzed  the  new  potential  of 
Nitz,  et  al .  (Ref.  38),  which,  as  mentioned  above,  is  similar  to  that  of  Foreman, 
et  al.  (Ref.  34). 

The  best  fit  of  our  calculated  intensities  to  the  experimental  data  of  the 
Hopfield  helium  continuum  requires  a  repulsive  energy  of  3.4  t  0.1  eV  for  the 
ground  state  of  Hej  in  the  vicinity  of  1.04  A  internuclear  separation.  This 
result  is  based  on  an  average  of  the  Foreman  and  Ceperley  potentials  since  the 
experimental  uncertainties  in  the  measured  intensity  distributions  prevent  us 
from  making  a  statistically  favored  choice. 
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The  results  of  this  study  now  cast  doubt  on  our  earlier  parameter izat ion 

(Ref.  23)  of  the  calculated  4JI  state  of  He«,  which  has  been  shown  to  be  the 

g  l 

upper  state  in  the  electron  autodetachment  experiments  of  Peterson,  et  al .  (Refs. 
20  and  21)  and  Kvale,  et  al.  (Ref.  22).  In  our  theoretical  study,  the  calculated 
location  of  the  metastable  4IIg  state  of  Hej  corresponded  to  a  predicted  electron 
affinity  of  0.182  eV,  in  good  agreement  with  the  experimental  estimate  of 
0.175  eV  reported  by  Kvale,  et  al.  (Ref.  22).  This  leaves  only  the  equilibrium 
internuclear  separation  of  He2  as  an  adjustable  parameter  to  reproduce  the 
experimental  vertical  separation  of  15.7  eV  between  the  lowest  vibrational  level 
of  the  4Ilg  state  of  Hej  and  the  ground  state  Hej  potential.  We  find  that  our 
predicted  equilibrium  separation  for  He^  (1.064  A)  must  be  shifted  outward  by 
~0.08  A  to  obtain  a  satisfactory  fit  to  the  autodetachment  data.  If  correct, 
this  represents  perhaps  the  first  known  instance  where  an  ab  initio  Cl  result 
predicted  too  small  a  value  for  the  equilibrium  separation.  Further  theoretical 
studies  of  He2  and  He^  ore  now  in  progress  using  a  complete  basis  set  (CBS) 
extrapolation  (Ref.  39). 


There  have  been  numerous  thermochemical  and  kinetics  studies  of  weakly  bound 
ionic  clusters  during  the  past  decade  (Refs.  40  and  41).  Of  particular  interest 
are  the  cation  and  H~  anion  sequences  since  they  represent  examples  of  the 
chemistry  of  strong  chemical  bonding,  in  H^,  the  simplest  hydrogen  bonded 
molecular  ion,  in  H^,  and  the  simplest  atomic  anion,  H“.  In  addition,  for  the 
first  few  members  of  these  sequences,  detailed  £b  initio  quantum  mechanical 
calculations  of  their  equilibrium  bond  lengths  and  vibrational  spectra  can  be 
carried  out . 

Very  recently,  the  possibility  of  forming  large  H*  clusters  (n  >  1000),  that 
could  be  kept  for  long  periods  of  time  in  a  Penning  trap,  has  been  discussed 
(Ref.  42).  Such  a  system  is  a  prototype  for  bulk  antiproton  or  antihydrogen 
storage,  with  uses  as  high  energy  density  material  or  as  an  energy  source  for 
advanced  propulsion  concepts.  The  mechanism  for  formation  and  storage  of  these 
large  cluster  ions  can  be  understood  only  if  the  thermodynamics  and  kinetics  of 
the  cluster  growth  processes  are  identified  and  analyzed  in  detail.  The  purpose 
of  our  study  was  to  examine  the  thermodynamic  stability  of  the  first  few  members 
of  the  hydrogen  cation  and  anion  cluster  sequences  since  their  stability  is  a 
requirement  for  defining  a  useful  association  pathway  for  the  growth  of  large 
clusters . 

Numerous  experiments  (Refs.  43  and  44)  have  measured  the  stability  of  the 
odd-merabered  hydrogen  clusters  that  are  formed  by  the  association  reaction: 


<  *  « 2  ♦ 


n  -  3,5,7, 


(10) 


10 


.'XN’.v.v, 
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These  clusters  exhibit  stabilization  energies  of  approximately  3-6  kcal/mol,  in 
good  agreement  with  theoretical  calculations  (Refs.  45  and  46).  Until  recently, 
the  even-membered  hydrogen  clusters  had  not  been  observed,  with  the  exception  of 
the  suggestion  of  a  stable  H£  ion  formed  by  the  attachment  of  a  H  atom  to  the 
ion  (Ref.  47).  A  very  recent  study  of  the  even-membered  clusters  (Ref.  48) 
suggests  that  the  H£  observation  could  be  interpreted  as  HjD*.  This  study  also 
indicates  a  high  degree  of  stability  for  H£. 

In  contrast,  the  negative  hydrogen  ion  clusters  have  never  been  observed  and 
a  recent  careful  search  for  stable  or  metastable  and  produced  negative 
results  (Ref.  49).  There  have  been  several  theoretical  studies  (Refs.  50  and  51) 
which  support  the  existence  of  stable  H3,  and  higher  clusters,  but  more  accurate 
ab  initio  calculations  (Ref.  52)  shed  doubt  on  this  conclusion. 

In  order  to  establish  the  thermodynamic  stability  of  these  hydrogen  cluster 
ions,  ab  initio  calculations  were  carried  out  to  determine  their  electronic 
structure  and  stable  geometries.  Several  Gaussian  basis  sets,  with  increasing 
flexibility,  were  employed  in  these  studies  and  levels  of  theory  up  to  single  and 
double  configuration  interaction  expansions  (CISD)  were  examined. 

The  results  of  these  calculations  are  summarized  in  Table  III,  which 
includes  the  H~  and  H2  constituents  as  well  as  H3.  It  is  clear  from  Table  III 
that  the  inclusion  of  the  vibrational  zero  point  energy  is  necessary  to  assess 
the  stability  of  H3.  We  find  with  every  basis  and  level  of  theory  that  the  H3 
ion  is  not  stable  with  respect  to  H~  and  H2 .  This  conclusion  is  supported  by  the 
results  of  earlier  theoretical  work  (Refs.  51  and  52),  and  by  the  recent 
experimental  study  of  Bae,  Coggiola,  and  Peterson  (Ref.  49). 

Our  calculations  appear  to  rule  out  the  formation  of  hydrogen  anion  clusters 
though  either  of  the  association  reactions: 


♦  H  *  Hn*  1 

(11) 

Hn 

♦  H2  *  HnZ2 

(12) 

Although  previous  studies  by  Hirao  and  Yamabe  (Ref.  52)  indicate  a  weak  stability 
for  H"  clusters,  up  to  n  *  13,  the  non-existence  of  and  H3  piecludes  formation 
of  higher  clusters  through  either  reaction  (11)  or  (12).  There  remains  the 
possibility  (albeit  remote)  that  association  of  H"  with  a  stable  Rydberg  state  of 
H  3  could  produce  H£.  Further  calculations  are  in  progress  to  access  this 
association  pathway. 
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PUBLICATIONS  AND  PRESENTATIONS 


The  significant  research  results  obtained  under  this  Contract  during  the 
past  year  have  been  prepared  for  publication  in  technical  journals  or  presented 
at  technical  meetings.  These  papers  and  meetings  are  listed  below.  Abstracts  of 
the  published  papers  are  included  in  the  Appendices  to  this  report. 

A.  Technical  Reports  in  Journals  and  Books 

1.  "Spectral  Intensity  of  the  Hopfield  Helium  Continuum:  An  Analysis  of  the 
Ground  State  Potential  for  He2,"  H.  H.  Michels,  R.  H.  Hobbs  and  J.  R. 
Peterson,  Chemical  Physics  Letters,  1 34 ,  1987,  pp.  571-574. 

2.  "The  Electronic  Structure  and  Stability  of  the  Anion,"  H.  H.  Michels  and 
J.  A.  Montgomery,  Jr.,  submitted  to  Chemical  Physics  Letter,  1987. 

3.  "Negative  Ion  Format  ion  in  Lithium  Atom  Collisions,"  H.  H.  Michels  and  J.  A. 
Montgomery,  Jr.,  to  be  published  in  the  Proceedings  of  the  Fourth 

Internet ional  Symposium  on  the  Production  and  Neutralization  o  f  Neg  ative  Ions 
and  Beams  (Brookhaven,  1986),  ed.  J.  Alessi,  1987. 

4.  "Electronic  Structure  and  Stability  of  Small  Cation  and  Anion  Hydrogen 
Clusters,"  J.  A.  Montgomery,  Jr.  and  H.  H.  Michels,  to  be  published  in  the 
Proceedings  of  the  1987  Cooling,  Condensation  and  Storage  of  Hydrogen  Cluster 
Ions  Workshop,  1987. 

B.  Technical  Paper  and  Lecture  Presentations 

1.  "Negat  ive  Ion  Formation  in  Lithium  Rydberg  Atom  Collisions."  Presented  at 
the  SDIO-Neutral  Particle  Beam  Annual  Program  Review,  Los  Alamos  National 
Laboratory,  New  Mexico,  September  30  -  October  2,  1986. 

2.  "The  Hopfield  Helium  Continuum  and  He^  (X  * £* )  Potential."  Presented  at  the 
Thirty-ninth  Annual  Gaseous  Electronics  Conference,  Madison,  Wisconsin, 
October  7-10,  1986. 

3.  "Assessment  of  the  He~  Potential  from  the  Autodetachment  Spectrum." 

Presented  at  the  Thirty-ninth  Annual  Gaseous  Electronics  Conference,  Madison, 
Wisconsin,  October  7-10,  1986. 

4.  "Neg  ative  Ion  Formation  in  Lithium  Atom  Collisions."  Presented  at  the  Fourth 
Int ernat ional  Symposium  on  the  Production  and  Neutralization  o  f  Neg  ative  Ions 
and  Beams,  Brookhaven  National  Laboratory,  Brookhaven,  Long  Island,  New  York, 
October  26-31,  1986. 
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5.  "Electronic  Structure  and  Stability  of  Small  Cation  and  Anion  Hydrogen 
Clusters."  Presented  at  the  1987  Workshop  on  the  Cooling,  Condensation  and 
Storage  of  Hydrogen  Cluster  Ions,  Menlo  Park,  California,  January  8-10,  1987 

6.  "Application  in  Computational  Chemistry:  Theoretical  Studies  of  Storable 
Highly  Energetic  Molecules."  Presented  at  the  Leermaker's  Symposium, 
Wesleyan  University,  Middletown,  Connecticut,  May  8,  1987. 

7.  "Assessment  of  the  Hej  Potential  from  the  Autodetachment  Spectrum." 

Presented  at  the  18th  Annual  Meeting  of  the  Division  of  Atomic,  Molecular, 
and  Optical  Physics,  Cambridge,  MA,  May  18-20,  1987. 
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TABLE  I 


LITHIUM  ATOM  INTERACTIONS 

(1)  Resonant  line  excitation 

Li l  2S ( 2s )  ]  +  hv  Li*  [2P(2p)  l(saturated)  ^2s-2p  *  671-°  001 

(2)  Rydberg  state  excitation 

Li* [ 2P ( 2p) ]  +  hv'  ♦  Li**[n>3]  Rydberg  state  excitation 

A2p-3s  *  812*8  nm 

(3)  Rydberg-ground  state  reaction 

Li[  2S(2s)  ]  +  Li**[n>3]  -►  Li^I 1 S]  +  Li'Ps] 

(4)  Rydberg-excited  2P  state  reaction 

Li*[ 2P(2p)  ]  +  Li**(n>3)  -►  Li+Ps]  +  Li"[lS] 

(5)  Rydberg-Rydberg  excited  state  reaction 

Li**(n>3)  ♦  Li**(n>3)  -  Li+^S]  ♦  Li^Sl 

jc’ic 

(6)  Associative  ionization  reaction  (loss  of  Li  ) 

Li l  2S (2s)]  +  Li**(n>4)  ♦  Li*  [X2!*]  ♦  e 

.  4rk  v 

(7)  Dissociative  ionization  reaction  (loss  of  Li  ) 

Lil2S(2s)]  +  Li**(n>4)  -►  Lit [A2Z*]  ♦  e  ♦  Li+[1S]  ♦  Lit2S]  +  e 
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TABLE  II.  Least-squares  vibrational  population  analysis  of  the  calculated 
Bopfield  spectral  distribution 
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Vibrational  Level  Ground  State  Re^  Potential 
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APPENDIX  A 


Spectral  Intensity  of  the  Hopfield  Helium  Continuum: 
An  Analysis  of  the  Ground  State  Potential  for  He^  . 

H.  H.  Michels  and  R.  H.  Hobbs 
United  Technologies  Research  Center 
East  Hartford,  CT  06108 


J.  R.  Peterson 
SRI  International 
Menlo  Park,  CA  94025 

ABSTRACT 


The  intensity  distribution  of  the  Hopfield  helium  continuum  has  been 
analyzed  using  accurate  numerical  wavefunct ions .  The  best  comparison  of  our 
calculated  intensity  to  the  experimental  continuum  of  Huffman,  et  al ,  yields  a 
repulsive  energy  of  3.4  i  0.1  eV  for  the  ground  X  state  of  He2  in  the 

vicinity  of  1.04  A  internuclear  separation. 


*Supported  in  part  by  AFOSR  under  Contracts  F49620-85-C-0095  and  F49620-85-R-0017 
and  NSF  Grant  PHY-84 10980. 

Published  in  Chemical  Physics  Letters,  134,  1987,  pp.  571-574. 
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APPFVDIX  B 


Negative  Ion  Formation  in  Lithium  Atom  Collisions 


H.  H.  Michels  and 
J.  A.  Montgomery,  Jr. 

United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 


ABSTRACT 


The  formation  of  Li-  by  dissociative  attachment  in  e  ♦  Li2  collisions  is 
characterized  by  a  large  cross  section  for  electron  attachment  to  highly  vibra- 
tionally  excited  Li2  molecules.  However,  the  electronic  structure  of  the  Li2~ 
system  dictates  that  low  energy  electrons  (<  0.2  eV)  are  needed  to  minimize  loss 
processes.  An  alternate  production  route  for  Li-  is  possible  via  direct  ion-pair 
formation  through  collisions  of  highly  excited  (Rydberg)  Li  atoms.  The 
mechanisms  for  ion-pair  formation  will  be  discussed  in  terms  of  the  high-lying 


molecular  Rydberg  states  of  Li 


2  ' 


♦Supported  in  part  by  AFOSR  under  Contract  F49620-85-C-0095. 

To  be  published  in  the  Proceedings  of  the  Fourth  International  Symposium  on  the 
Production  and  Neutralization  of  Negative  Ions  and  Beams  (Brookhaven,  1986),  ed . 
J.  Alessi,  1987. 
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APPENDIX  C 

Electronic  Structure  and  Stability  of  Small  Cation 
and  Anion  Hydrogen  Clusters 


J.  A.  Montgomery,  Jr.  and  H.  H.  Michels 
United  Technologies  Research  Center 
East  Hartford,  CT  06108 


ABSTRACT 


Ab  initio  calculations  of  the  electronic  structure  of  H  ”  and  H  *  clus- 
-  n  n 

ters  have  been  carried  out  using  accurate  Gaussian  basis  sets  and  with  levels 
of  theory  up  to  fourth-order  perturbation  theory  (MPA)  and  single  and  double 
excitation  configuration  interaction  (CISD).  The  odd  hydrogen  cation  addition 
sequence,  Hn+  ♦  H2  ♦  ^n+2*  aPPears  to  thermodynamically  stable  for  large 
size  cluster  formation.  The  even  H&+  cation  also  exhibits  surprising  stabil¬ 
ity  in  Djj  symmetry.  In  contrast,  the  hydrogen  anion  addition  sequences,  Hn" 

♦  H  ♦  Hn"j  and  Hn"  ♦  H2  ♦  Hn"2»  appear  to  be  thermoneutral  or  unstable. 


♦Supported  in  part  by  AFOSR  under  Contract  F49620-85-C-0095  and  in  part  by 
AFRPL  under  Contract  F04611-86-C-0071 . 

To  be  published  in  the  Proceedings  of  the  1987  Cooling,  Condensation  and 
Storage  of  Hydrogen  Cluster  Ions  Workshop,  1987. 
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APPENDIX  D 

The  Electronic  Structure  and  Stability  of  the  Hj  Anion* 

H.  H.  Michels  and  J.  A.  Montgomery,  Jr. 

United  Technologies  Research  Center,  East  Hartford,  CT  06108,  USA 


ABSTRACT 


A  systematic  study  of  the  electronic  structure  of  the  anion  has  been 
carried  out  using  both  perturbation  theory  and  configuration  interaction. 
Counterpoise  corrections  are  used  to  estimate  the  basis  set  superposition  error 
in  the  computed  dissociation  energy.  We  find  that  is  not  thermodynamically 
stable  relative  to  H”  ♦  H2  upon  consideration  of  the  vibrational  zero-point 
contribution  to  the  dissociation  energy.  Isotopic  analysis  predicts  weak 
stability  for  and  Dj-H-. 


$ 


v. 


♦Supported  in  part  by  AFOSR  under  Contract  F49620-85-C-0095 .  To  be  published  in 
Chemical  Physics  Letters,  1987. 
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